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\ We present a new set of kinematic equations for front motion in bistable media. 

The equations extend earlier kinematic approaches by coupling the front curvature 
with the order parameter associated with a parity breaking front bifurcation. In 
addition to naturally describing the core region of rotating spiral waves the equa- 
tions can be be used to study the nucleation of spiral-wave pairs along uniformly 
propagating fronts. The analysis of spiral- wave nucleation reduces to the simpler 
problem of droplet, or domain, nucleation in one space dimension. 



1 Introduction 



The onset of spatio-temporal disorder in reaction-diffusion systems is often 
accompanied by the spontaneous nucleation of spiral waves. The conditions 
^ ■ for spiral- wave nucleation, mechanisms behind it, and implications on the 

resulting dynamics have been studied extensively [1-6]. Yet detailed studies 
of the nucleation process itself are still lacking. One difficulty in carrying out 
such studies stems from the two-dimensional structures of spiral waves. In this 
paper we report on the derivation of new kinematic equations for front motion 
in bistable media which reduce the two-dimensional spiral-wave nucleation 
problem to a one-dimensional droplet nucleation problem [7,8]. We use the 
kinematic equations to demonstrate the destabilization of traveling V-shape 
solutions to spiral nucleation. 

The derivation is carried out in a parameter range that includes a pitchfork 
front bifurcation called a nonequilibrium Ising-Bloch (NIB) bifurcation [9- 
12]. At a NIB bifurcation, a stationary "Ising" front loses stability to a pair 



1 aric@lanl.gov 

2 ehud@bgumail.bgu.ac.il 



Preprint submitted to Elsevier Preprint 



9 February 2008 




Fig. 1. (a) The nonequilibrium Ising-Bloch (NIB) front bifurcation. The solid line 
represents a branch of stable front solutions with speed c. At the bifurcation point, 
a = a c , the stationary solution becomes unstable to a pair of counterpropagating 
fronts (b) The NIB front bifurcation and planar front transverse instability bound- 
aries in the e — 5 parameter plane. The thick line is the NIB bifurcation, <5f(c), 
and the dashed lines are the boundaries for the transverse instability of Ising, Si(e), 
and Bloch, 5s (e), fronts. Above these lines planar fronts are unstable to transverse 
perturbations. The thin lines are the approximations to the transverse instability 
boundaries obtained from the kinematic equations. Parameters: a\ = 4.0, ao = 0. 

of counter-propagating "Bloch" fronts as a parameter a is decreased past a 
critical value (Fig. la). The coexistence of the two Bloch fronts beyond the 
bifurcation allows for nucleation events where a front line segment of one Bloch 
type undergoes a transition to the other Bloch front. 

The kinematic description of front motion consists of the following equations: 

• An equation for the order parameter, Co, associated with the NIB bifurca- 
tion: 

^ = (a c -a)Co-f3C* + 1 K + lo + ^-^J ' nC n ds' . (1) 

o 

• A geometric equation for the front curvature, k: 



• An equation relating the normal front velocity C n , the curvature k, and the 
order parameter Co: 

C n = C -Dk. (3) 

In these equations s is the front arc length and the critical parameter value 
a c designates the NIB bifurcation point. The kinematic equations are derived 
for a FitzHugh-Nagumo model with a diffusing inhibitor. They generalize an 



2 



earlier kinematic approach [13,14] by treating Co as an independent dynamic 
mode rather than a constant. 



The order parameter equation (1) yields the bifurcation diagram of Fig. la for 
planar (not curved) front solutions in a symmetric system (70 = 0). The two 



Bloch branches pertain to high and low C values, C = C$ = ±y (a c — a)/ (3. 
Because of the coupling between the order parameter and curvature equations, 
curvature variations are capable of nucleating low Co segments (droplets) in 
regions of high Co- Droplet nucleation in the kinematic equations corresponds 
to spiral-wave pair nucleation in the physical two-dimensional plane. 

Spatially uniform front transitions induced by curvature, producing for exam- 
ple breathing spots, have been studied in Ref. [15,16]. The present paper is an 
extension of this earlier work to nonuniformly curved fronts. Complementary 
results to those reported here appear in Ref. [17]. 



2 Derivation of the kinematic equations 

The curvature equation (2) follows from purely geometric considerations [13,14,18- 
20]. The normal velocity relation (3) and the order parameter equation (1) can 
be derived using the FitzHugh-Nagumo reaction-diffusion model with a dif- 
fusing inhibitor, 



where u and v, the activator and the inhibitor, are real scalar fields and V 2 
is the Laplacian operator in two dimensions. The parameter a\ is chosen so 
that (4) describes a medium with two stable spatially uniform states: an up 
state (u+,v+) and a down state (w_,t>_). Ising and Bloch front solutions con- 
nect the two uniform states (u±, v±) as the spatial coordinate normal to the 
front goes from —00 to +00. The remaining parameter space is spanned by 
e, 5 and do, or alternatively by rj = y/e5, \x = e/5, and a . Note the parity 
symmetry (u,v) — > (—u, —v) of (4) for a = 0. 

The NIB bifurcation line in the e — 5 plane for ao = is shown in Fig. lb. 
For fi <C 1 it is given by 5 = 5f(c) = Vc/ e > or V = Vc, where i] c = 2 ^ q3 and 
q 2 = ax + 1/2 [11]. The single stationary Ising front that exists for 77 > r/ c loses 
stability to a pair of counter-propagating Bloch fronts at rj = i] c . Beyond the 
bifurcation (rj < i] c ) a Bloch front of an up state invading a down state coexists 
with another Bloch front of a down state invading an up state. Also shown in 
Fig. lb are the transverse instability boundaries (for a = 0), 5 = 5i(e) = e/r/ 2 




du 

dt 
dv 

dt 



u — a\v — ao + V 2 v , 



e^(u - u 3 - v) + 5 _1 V 2 m, 



(4) 
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and S = ^s(e) = rj c /-^e, for Ising and Bloch fronts respectively. Above these 
lines, 5 > Sj t B, planar fronts are unstable to transverse perturbations [6,21]. 
All three lines meet at a codimension 3 point P3: e — rj 2 , 5 — 1, ao = 0. 

The following assumptions are used in deriving Eqs. (3) and (1): /i = e/8 <C 1, 
the front speed c is small, the front is weakly curved, and curvature variations 
along the front are weak. The second assumption is met by considering nearly 
symmetric (|a | << 1) systems and restricting parameter values to the Ising 
regime or the vicinity of the NIB bifurcation where the front speed is c oc 
a/?7 c — i] <C 1. The last two assumptions are met by considering parameter 
values below or just beyond the threshold for the transverse instability of 
fronts. 

The kinematic equations are derived by transforming to an orthogonal coordi- 
nate system that moves with the front and using singular perturbation theory 
with /i = e/5< 1 asa small parameter [15]. The analysis of the narrow front 
core region, where u changes on a scale of order y/Jl, gives equation (3) with 

Co = %=v f (s,t), D = 5~\ (5) 



where Vf(s,t) is the approximately constant value (in the direction normal to 
the front) of the inhibitor v in this region. 

Away from the front core, where both u and v vary on a scale of order unity, a 
free boundary problem is obtained by using the solutions u±(v) of u—u 3 —v = 
in the inhibitor (v) equation. Matching the solutions from each side of the front 
at the front core (in the limit /i — > 0) gives the equation 



dv f _ V2(r] c - 7]) 



dt 



QTc 
2(1 



v f 



3q 



-K + 



4?7, 

d 2 Vf 
ds 2 



3°° 
ds 



(6) 



Equation (6) coincides with (1) using the following identifications: C = 

~7^2 V f> a = if' ac = £' P = 1 / 6 ' 7 = ac ^ ~~ 6 ~^> and 70 = 2a ^- 
More details about the derivation appear in Ref. [17]. 



3 Numerical solutions of the kinematic equations 
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Fig. 2. A front solution to the kinematic equations (l)-(3). (a) The order parameter 
Co and the curvature k along the arc length s. (b) In the x—y plane the front solution 
corresponds to a rotating spiral wave. The shaded (light) region corresponds to an 
up (down) state. Parameters: a\ = 4.0, clq = 0, e = 0.01234, S = 1.0. 

3.1 Spiral waves 

Consider a front solution in the kinematic equations that connects the state 
Cq = Cq, k = 0, at s = — oo with the state Co = Cq , k — 0, at s — +oo, where 
for a symmetric model (a = or 7 = 0) C$ = ±^J (a c — a)/ (3. Fig. 2a shows 
such a solution obtained by numerically integrating (l)-(3). As demonstrated 
in Fig. 2b this front solution of the kinematic equations (l)-(3) represents 
a spiral-wave in the FitzHugh-Nagumo model (4). Unlike earlier kinematic 
approaches [13,14] the spiral core is naturally captured. 

3.2 Traveling V-shape solutions and spiral-wave nucleation 

For Cq constant, the curvature equation (2) and the normal velocity rela- 
tion (3) have a family of traveling V-shape solutions [20,22] 

<s) = - — *f b 2 = C n (0) 2 - Cl , (7) 

Co + C n (0) cosh(os) 

where C n (0) = C — Dk(0) is an arbitrary constant. The traveling V-shape 
solution (7) is also an exact solution of the kinematic equations (l)-(3) for 
5 = 1. In that case 7 = a c (l — 5 _1 ) = and the order parameter equation (1) 
decouples from the curvature equation for constant Co solutions. The specific 
values C can assume are determined as solutions of the cubic equation (a c — 
a)C - PCI + 70 = 0. 
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Fig. 3. Nucleation of a spiral-wave pair from an initial traveling V-shape front. Left 
column: the Cq(s) and k(s) profiles. Right column: the front line shape in the x — y 
plane. Parameters: a x = 4.0, a = -0.0001, e = 0.0115, S = 1.063. 

We have also found stable traveling V-shape numerical solutions for 5 ^ 1. 
Fig. 3a shows a V shape Bloch front traveling stably at constant speed. When 
approaching the NIB bifurcation (by increasing e and/or 5) this solution be- 
comes unstable and nucleates spiral-wave pair. The time evolution of an un- 
stable V-shape solution is shown in Fig. 3(b-d). 

We carried out numerical solutions of the FitzHugh-Nagumo model (4) to 
test the prediction, based on solutions of the kinematic equations, that trav- 
eling V-shape solutions are destabilized as the NIB bifurcation is approached. 
Choosing parameter values close to those of Fig. 3 an initial V form is indeed 
unstable and a pair of spiral- waves nucleate in the region of highest curvature 
as shown in Fig. 4. After the spirals form the resulting traveling fronts may 
either repel each other or reconnect and the V-shape solution will continue 
to propagate with the possible nucleation of new spiral pairs. Interactions 
between front segments are not included in the kinematic equations (l)-(3) 
so the resulting evolution of the spiral pair and possible front reconnections 
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Fig. 4. Spiral-wave nucleation from an unstable traveling V-shape solution of the 
FitzHugh-Nagumo equations (4). The initial traveling V-shape solution (a) slows 
down (b) at the high curvature region and nucleates a pair of spiral waves (c). 
Depending on the parameters the resulting fronts may either repel or reconnect (d). 
Parameters: e = 0.0122, S = 1.064, a x = 4.0, a = -0.0001. 

cannot be captured. Front interactions have been studied recently in the fast 
inhibitor limit (rj ^> i] c ) [23]. 



4 Conclusion 

The kinematic equations presented in this paper extend an earlier kinematic 
approach using the curvature equation (2) and a linear relation between the 
normal velocity and curvature [13,14]. Indeed, far into the Bloch regime (a 
significantly smaller than a c ), Co is no longer a slow degree of freedom and can 
be eliminated adiabatically. Equations (1) and (3) then reduce to a multivalued 
algebraic relation Co = Co(/t) whose three branches (two Bloch and one Ising) 
give approximately linear C n - k relations [21]. 

But close to the NIB bifurcation C becomes an active degree of freedom and 
the earlier approach, often called the "geometric approach" , is no longer valid. 
The new kinematic equations capture the core structures of spiral waves and 
spiral-wave nucleation processes. Since one space dimension has been elim- 
inated in deriving the kinematic equations, the two-dimensional problem of 
spiral-wave nucleation has been reduced to a simpler droplet nucleation prob- 
lem in one space dimension. 
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